Effect of running coupling on photon emission from quark gluon plasma 
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We discuss the role of running coupling on the thermal photon yield from quark gluon plasma. 
It is shown that the photon production rate from the partonic phase is considerably enhanced 
when running coupling is considered with respect to a fixed value. However, we show by explicit 
evaluation that although this difference survives the space-time evolution the experimental data 
cannot distinguish between the two once the hard contribution, which is an essential component of 
photon production mechanism, is added. 



Detection of quark gluon plasma (QGP) in heavy 
ion collisions has received significant attention in recent 
years. Among possible signals, electromagnetic probes 
arc one of the most promisingtools to characterize the 
initial state of the collisions [lj. Owing to their weak 
coupling with the constituents of the system they tend 
to escape almost unscattered. In fact, photons (dilcptons 
as well) can be used to determine the initial temperature, 
or equivalently the equilibration time. By comparing the 
initial temperature with the transition temperature from 
lattice QCD, one can infer whether QGP is produced. 

Photons are produced at various stages from (i) initial 
hard scattering of partons, (ii) scattering of charged par- 
ticles in the thermal medium (QGP and hadronic matter) 
and (iii) from ir° and 77 decays. If this decay contribu- 
tion is subtracted from the total photon yield what is left 
is the direct (excess) photons. The thermal photon rate 
due to Compton and annihilation processes in a quark 
gluon plasma has been calculated by several authors over 
the last two decades In all these calculations the 

strong coupling, a s is treated as constant or function of 
temperature T. However, in case of rclativistic heavy ion 
collisions, apart from temperature there is also the mo- 
mentum scale k. One has to take into account the case 
when k ~ T and treat a s to be function of both k and 
T • By incorporating this fact it is shown that the en- 
ergy loss is a factor of 2 — 6 more than the case when a s 
is constant @, HI ■ The energy loss calculation using run- 
ning coupling and reduced screening mass in Q , explains 
single electron Raa quite well. It is the purpose of this 
brief report to treat the strong coupling as running and 
apply it to the case of thermal photon production from 
QGP. 

The lowest order processes for photon emission from 
QGP are the Compton scattering (q(q) g —> q(q) 7) and 
annihilation (qq —> 177) process. The total cross-section 
diverges in the limit t or u — > 0. These singularities have 
to be shielded by thermal effects in order to obtain infra- 
red safe calculations. It has been argued in Ref. @ that 
the intermediate quark acquires a thermal mass in the 
medium, whereas the hard thermal loop (HTL) approach 
of Ref. 0] shows that very soft modes are suppressed in a 
medium providing a natural cut-off k c ~ gT. We assume 



that the singularities can be shielded by the introduc- 
tion of thermal masses for the participating partons. The 
differential cross-sections for Compton and annihilation 
processes are given by (loj . 
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where m is the in-medium thermal quark mass, m 2 = 
ra t h 2 = 27ra s T 2 /3, a e and a s are the electromagnetic 
fine-structure and the strong coupling constants respec- 
tively. The static photon rate in 1 + 2— s-3 + 7 can be 
written as [l| 
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.Mi represents the amplitude for Compton or annihila- 
tion process. The overall degeneracy factor Ni ~ 320/3 
and 20 for Compton and annihilation processes respec- 
tively involving u and d quarks. 

As mentioned earlier, the infra-red cut-off is fixed by 
plasma effects, where only the medium part is considered, 
completely neglecting the vacuum contribution leading 
to ambiguity in the calculation of cross-section at finite 
temperature QCD. If the latter part is taken into ac- 
count the strong coupling should be running. Thus for 
any consistent calculation one has to take this fact into 
consideration. We have in that case a s = a s (k, T) where 

k = VW\- 

Photons from thermal hadronic matter also make up 
an essential component of the total photon yield from 
heavy ion collisions. These are emitted in reactions be- 
tween charged hadrons and in the radiative decays of 
unstable hadrons [l], [TJ Ell • I n this work, we have used 
the amplitudes of photon producing reactions involving 
the 7r, p, w, 77, K and K* mesons obtained in (T3 |. 

The hard photon contribution can be calculated by 
perturbative QCD (pQCD). In order to calculate pho- 
ton production from reactions of the type \ia hs — > 7 X 
(where h,A,h,B refer to hadrons), we assume that the en- 
ergy is such that the partonic degrees of freedom be- 
come relevant and they behave incoherently. The cross- 
section for this process can then be written in terms of 
elementary parton-parton cross-section multiplied by the 
partonic flux which depends on the parton distribution 
functions [l4[. The energy scale for this to happen i.e. 
the factorization scale is denoted by Q 2 , the square of 
the momentum transfer of the reaction. Starting with 
two body scattering at the partonic level the differential 
cross-section for the reaction of above type can be written 
as Gl 
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where, xt = 2pr/ys and the factor K is introduced to 
take into account the higher order effects. A few com- 
ments about the K factor is in order here. The cross- 
section in the above expression is calculated perturba- 
tively to leading order (LO) in the strong coupling. In 
cases where the next-to-leading (NLO) order terms are 
comparable to the LO terms the K factor defined as 
NLO/LO is introduced in the LO computations to bring 
in the essence of the NLO terms. It has been shown 
in [l6j that K depends on the choice of the momentum 
scale, the parton distribution functions and the shadow- 
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FIG. 1: Strong coupling as a function of momentum scale at 
two different temperatures Q. 



ing effect and its value lies between 2-3. In the present 
calculation we take K ~ 2.5. 

We also include photons from fragmentation process. 
This is accomplished by introducing the fragmentation 
function, £> 7 / c (z, Q 2 ), which when multiplied by dz gives 
the probability for obtaining a photon from parton c, z 
being the fractional momentum carried by the photon. 
Once the photon production cross-section is obtained 
from hadron-hadron collision we can now determine the 
direct photon production rates due to hard scattering 
between partons from nucleus- nucleus collisions at rela- 
tivistic energies. To do this we must note that the ex- 
perimental data are given for a particular ccntrality. In 
order to take this into account we introduce the central- 
ity parameter which depends on the maximum impact 
parameter b m . The photon yield from hard collisions is 
then calculated from the expression: 
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where TZ(b m ) = (ABTab) = -Fr m 

and TAB(b) = / gPsTa(s) Te(b — s) is the nuclear over- 
lap function. For - 10 % centrality we obtain 1Z ~ 
21.7m6 -1 . Before going to the numerical evaluation of 
the static photon rate we plot the running coupling in 
Fig. ([1]) where the parametrization for a s (k,T) is taken 
from Ref. Q. It is seen that the value of the coupling is 
largest when k ~ T. For k >> T it agrees well with the 
temperature dependent coupling. These features of a s 
have important consequence on the photon production 
rate as we shall see below. 

The static photon rate is obtained from Eq.(3) using 
the running coupling. For T = 200 MeV the rates are 
shown in Fig. ([2]). The photon emission rate is enhanced 
by a factor of 1.7 - 6 compared to the case when the mo- 
mentum dependence of the strong coupling is neglected. 
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FIG. 2: Static photon rate with and without the running 
coupling. 



Photons are produced at all stages of the collision and 
so it is necessary to integrate the emission rates over the 
space-time volume from creation to freeze-out. We as- 
sume that quark gluon plasma having a temperature Tj 
is produced at an initial time Tj. Hydrodynamic expan- 
sion and cooling follows up to a temperature T c where 
QGP crosses over to a hadronic gas. Subsequent cool- 
ing leads to freeze-out of the fluid element into observ- 
able hadrons. In the present work the fireball is taken 
to undergo an azimuthally symmetric transverse expan- 
sion along with a boost invariant longitudinal expansion. 
The local temperature of the fluid element and the as- 
sociated flow velocity as a function of the radial coordi- 
nate and proper time is obtained by solving the the en- 
ergy momentum conservation equation T^ v = where 



g^ v p jg t ne energy momentum ten- 



sor for ideal fluid. This set of equations are closed with 
the Equation Of State (EoS); typically a functional rela- 
tion between the pressure P and the energy density e. It 
is a crucial input which essentially controls the profile of 
expansion of the fireball. To minimize model dependen- 
cies we take the EoS from the lattice calculations of the 
Wuppertal-Budapest Collaboration [l7j . 

The initial temperature is related to the experimen- 
tally measured hadron multiplicity through entropy con- 

servation as T?(b m )n = 45C(3) 2 ;^4 nfc (ffiftO) 
where (dN/dy(b m )) is the hadron (predominantly pions) 
multiplicity for a given centrality class with maximum 
impact parameter b mi Ra is the transverse dimension 
of the system and afc is the degeneracy of the system 
created. The hadron multiplicity resulting from Au + Au 
collisions is related to that from pp collision at a given im- 
pact parameter and collision energy through the relation 
(%(b m )) = [(1 - x){N part {b m ))/2 + x(N coll (b m ))} 
where x is the fraction of hard collisions. (N part ) is the 
average number of participants and {N co u} is the average 
number of collisions evaluated by using Glauber model. 
dN^/dy = 2.5 - 0.25 Ins + 0.023 In 2 s is the multiplic- 
ity of the produced hadrons in pp collisions at centre 
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FIG. 3: Thermal photon p T distributions for Ti = 400 MeV 
and Ti — 0.2 fm/c at central rapidity. Upper panel corre- 
sponds to vq — 0. Individual contribution from hard and 
fragmentation photons has been shown in the lower panel for 
vq = 0.2. 



of mass energy, y/s [l9j . We have assumed that 20% 
hard (i.e. x = 0.20 ) and 80% soft collisions are re- 
sponsible for initial entropy production. Considering 
- 10 % centrality we get % = 400 MeV for Ti = 0.2 
fm/c. One also requires the initial energy density and ra- 
dial velocity profiles which are taken respectively as (20j 



and v{ri,r) = vq 



1 



l +e (r-H A )/« 



1 + e (r-R A )/5 

where S (~ 0.5 fm) is a parameter, known as the surface 
thickness. As discussed in [2(J, this choice of the initial 
fluid velocity profile is motivated by the fact that for a 
physical system the initial fluid velocity is zero inside the 
matter which approaches a value vq which is of the or- 
der of a typical particle transverse velocity in the diffuse 
region outside the matter distribution. 

The other inputs are the transition temperature T c 
which is taken as 175 MeV as obtained from lattice 
QCD [2l], [22| and the freeze-out temperature, Tf which 
is taken to be 120 MeV. We now plot the thermal photon 
yield from both QGP and hot hadronic matter in Fig. ([3]) 
for vq — and 0.2 in the upper and lower panels. At 
very low pt ~ 0.5 GeV, the contribution from hadronic 
matter dominates. Beyond that the QGP contribution 
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FIG. 4: Photon p t distributions for Ti = 400 MeV and n = 
0.2 fm/c at central rapidity. The data (for \y\ < 0.35) are 
taken from Ref. [23| . 

starts to take over. Interestingly, the effect of running 
coupling on the thermal photon production from QGP 
does survive the space time evolution and continues to 
be discernible in this pt domain. The effect of non-zero 
initial velocity (vq) is also visible in the upward shift of 
the spectra at higher px in the lower panel compared to 
the upper. Moreover, the relative separation between the 
QGP contributions with and without the running cou- 
pling appears to be independent of the space-time evolu- 
tion scenarios corresponding to vq = and 0.2. To assess 
their relative importance in comparison to the thermal 
yield we also show in the lower panel the contribution 
from hard QCD photons and photons produced from jet 
fragmentation. As shown by the dash-dotted line this 
contribution clearly dominates the photon yield beyond 
about 2 GeV. Below that the hard and fragmentation 
contributions arc almost similar to QGP contribution. 



We now compare the total yield with the direct pho- 
ton data from Au+Au collisions at RHIC obtained by 
PHENIX in fig. (gj). It is observed that the data 
is best reproduced by assuming a small initial veloc- 
ity of the order of 0.2 (solid line) compared to vq = 
(dashed line). However, the curves for a s = o: s (T) and 
a s = a s (k, T) have merged with each other implying that 
the observed difference seen in the thermal photon contri- 
bution for the two cases has been washed away once the 
hard and fragmentation contributions arc added. Such 
a result can be understood once we realize that the con- 
tribution at a given value of the transverse momentum, 
especially up to 2-3 GeV/c is a superposition of contribu- 
tions from QGP at temperatures from Ti to T c , hadronic 
matter from T c to Tf as well as from hard scatterings. 
Although the QGP contribution clearly dominates for 
Pt > 0.5 GeV observation of momentum-dependence of 
the strong coupling in the transverse momentum spectra 
of single photons does not appear to be feasible as it is 
overshadowed by the contributions coming from initial 
hard collisions. 

To summarize, we have calculated the static photon 
rate from QGP due to Compton and annihilation pro- 
cesses using the temperature and momentum dependent 
strong coupling. The rate is then contrasted with the 
case where a s depends only on the temperature of the 
system. It is found that the static photon rate enhances 
significantly if the running coupling is used. We then 
perform a space-time evolution using rclativistic hydro- 
dynamics with initial conditions for Au+Au collisions at 
200 GeV/n at RHIC. The significant difference in the 
yields does survive the space time evolution and can be 
observed in the thermal photon spectra. However, due 
to the large contributions coming from the initial hard 
collisions and jet fragmentation the single photon data 
from PHENIX cannot distinguish between the scenarios 
with and without momentum dependence of the running 
coupling. 
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